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1 . INTRODUCTION 

In electronic equipment in which high density mounting 
is proceeding, a thermal fatigue crack problem of solder joint 
becomes more important. Thus, many studies have been made as 
to reliability evaluation of solder joint (1) " (4) . The authors 
also have proposed an elastic-creep model in which attention 
is paid to non-elastic behavior of Sn63-Pb37 solder and have 
indicated the effectiveness of the model {5> " {10) . Many of the 
evaluation methods, which have been proposed, are based on the 
concept of crack initiation lifetime determined by a 
Cof f in-Manson type fatigue lifetime prediction formula. An 
increase in electric resistance value, which is a problem in 
reliability, generally occurs at a time near to breaking 
lifetime. Accordingly, prediction of crack initiation 
lifetime is safety side evaluation, which is reasonable as a 
design method. However, knowledge as to crack evolution may 
be necessary under various design conditions. From the above 
background, studies on crack evolution also have been recently 
increased " (13 > . 

The authors also have reported an experimental method, 
in which a simple method (8) based on material damage law is 
combined with a fatigue test of actual joint (10), in addition 
to evaluation of crack evolution (13) based on breakdown 
dynamics. This article will introduce first a result of 
mechanical fatigue test of solder bump joint used in Ball Grid 
Array package (BGA) and the like, prediction of crack initiation 
lifetime, and a simple prediction method of breakdown lifetime . 
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In addition to the above, this article will reports creep J 
integration range determined by FEM analysis using the 
elastic-creep model (5) and a result of investigations as to crack 
evolution characteristics in bump joint mode. 

2. MECHANICAL FATIGUE TEST OF SOLDER BUMP JOINT AND LIFETIME 
PREDICTION 

2 . 1 Test method 

Fig. 1 shows the outline of a tester and an Sn63-Pb37 
solder bump joint specimen. In the test, the fatigue tester, 
which was driven by a linear actuator and a reduced displacement 
mechanism was developed and used. The specimen had a shape 
approximately the same as that of a BGA joint. The joint was 
subjected to shear deformation by repeatedly applying relative 
replacement to ceramic substrates between which the bump joint 
was sandwiched. The test was carried out using displacement- 
controlled symmetrical triangular waves at room temperature. 
Since the solder joint of an actual semiconductor package had 
a low strain rate, a low wave frequency had to be set in the 
fatigue test. This article will report a result of test carried 
out under the conditions of total displacement amplitude AX = 
7 to 14 ^in and wave frequency f «3 to 5xl0" 3 Hz. 

Fig. 1 Outline of mechanical fatigue tester and solder bump 
joint specimen 

2.2 Prediction of crack initiation lifetime 
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Fatigue cracks initiated approximately at the same time 
from the overall peripheries of the solders in the vicinities 
of the joint boundaries of the bump specimen, and the initial 
cracks developed into final breakdown. In the meantime, a 
maximum load range was monotonously dropped as the number of 
cycles increased. Fig. 2 shows an example of reduction of a 
load range. 

Fig. 2 Drop of load range 

In the solder bump fatigue test, fatigue crack cannot be 
observed from the outside. The authors propose to indirectly 
prescribe the number of cycles in which initial crack occurs 
from the reduction rate of the load range. Drop of the load 
range in the fatigue test is caused by drop of rigidity caused 
by initiation and evolution of crack and to the repeated 
softening phenomenon of Sn63Pb37 solder material. First, the 
drop of the load range caused by the drop of rigidity was 
determined using an FEM model having crack. As a result, the 
load range was dropped by 11% with respect to that of an initial 
shape model. Next, when the drop of the load range caused by 
the repeated softening was examined by a result of fatigue test 
of bulk member and an average shear strain range caused in bump 
portion, it was calculated that a load range drop rate was 7%. 
A load drop rate at the time of initiation of crack was obtained 
by multiplying the load ranges, which were calculated at the 
respective times at which crack initiated and softening was 
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repeatedly carried out after load was dropped. In the bump 
shape in the test, it was calculated that the load range drop 
rate when crack initiated at the depth of 50 jam was approximately 
constant and 17%. The number cycles in which fatigue crack 
initiated was determined based on the above calculation and 
arranged as shown in Fig. 3. 

FEM analysis of the bump joint specimen was carried out 
using the general-purpose structural analysis code ABAQUS based 
on a finite element method . The elastic-creep model {5) proposed 
by the authors was used as the characteristics of solder 
material. Fig. 4 shows a contour view of equivalent creep 
ranges. In the bump joint mode, since strain-concentrated- 
portion is located at a singularity point, lifetime prediction 
based on the maximum value of an angular point is improper. The 
authors proposed previously a method of providing an evaluation 
point within 50 jam from the singularity point. The authors 
calculated predicted lifetime by the following Cof f in-Manson 
equation (7> by applying the method. 

N f =0.146 • Ae m A 94 (1) 
where, N f shows the number of lifetime cycles, and Ae in shows 
a non-elastic strain range (here, equivalent creep strain 
range) . Fig. 3 also shows a result of lifetime prediction. 
Since bump specimens cannot help being dispersed in the shape 
thereof, crack lifetimes are different several times. However, 
this tendency excellently corresponds to the result of lifetime 
prediction based on the equivalent creep strain range obtained 
50 jam inside of the singularity point. 
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Fig. 3 Crack initiation lifetime 

Fig. 4 Distribution of equivalent creep strain range 

2.3 Simple prediction method of breakdown life 

The number of cycles in which the load range drop rate 
reached 75% was set as a breakdown lifetime. Fig. 5 shows the 
relationship between the total displacement amplitude AX and 
the number of cycles in which the load range drop rate reached 
17% and 75%. Although it cannot be said that the distribution 
of fatigue crack initiation lifetime (load drop rate: 17%) is 
in parallel with that of breakdown lifetime (load drop rate: 
75%) , both of them incline comparatively similarly. In the bump 
shape subjected to the test, it is clarified that breakdown 
initiation lifetime is about four to six times longer than crack 
initiation lifetime. A method of multiplying predicted crack 
initiation lifetime by the above ratio is contemplated as a 
simple prediction method of bump breakdown lifetime in actual 
design. 

Fig. 5 Relation between crack initiation lifetime and 
breakage lifetime 

3. ANALYSIS OF CRACK DEVELOPMENT 

3.1 Calculation of crack evolution rate 

Calculation similar to that described in the previous 
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paragraph was carried out by changing the depth of initial crack 
by several micron meters, and the relationship between crack 
depth and load range drop rate was examined. As a result, 
relationship between crack evolution rate and load range 
reduction when crack initiated at the depth of 50 was 2.7 
jam/%. Fig. 6 arranges the evolution rate immediately after the 
initiation of crack based on the above result. To verify the 
reasonability of the above prescribing method, breakdown was 
checked after the test was stopped in the predetermined number 
of cycles, and the state of crack evolution was observed. As 
a result, about 50 jun deep initial cracks occurred in the 
periphery of a bump under the condition of AX = 12 jam after 150 
cycles passed. It was observed that the cracks developed to 
a depth of about 100 pm after 300 cycles passed. An average 
evolution rate was 0 . 3 jam/cycle during the about period of time, 
which is in good agreement with the result arranged in Fig. 6. 

Fig. 6 Clack evolution rate immediately after occurrence of 
initial crack 

3.2 Calculation of AJc by elastic-creep analysis 

The FEM analysis of bump joint specimen was carried out 
using the general-purpose structural analysis code ABAQUS based 
on the finite element method. A three-dimensional analysis 
model was created as shown in Fig. 7 in consideration of the 
symmetrical property of the bump specimen. Cracks having a 
depth of 50 \im was formed around the entire periphery of the 
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vicinity of a joint interface. The elastic-creep model was used 
as the solder material characteristics. The condition of the 
test resides in a strain rate region (about 3 x 10" 4 /sec) to 
which the elastic-creep model can be applied. Similarly to the 
test, displacement-controlled triangular waves were applied to 
the upper surface of model as a load condition. A creep J 
integrated value was determined at point A in Fig. 7 in the 
direction of an arrow. For the convenience of an amount of 
calculation, the integrated value was calculated up to a fourth 
path. Fig. 8 shows an example of J* in a half cycle on a crack 
opening side. Since the crack opening side is in a path in the 
vicinity of a crack extreme end, a steady state is reached in 
a short time. A creep J integration range AJc was calculated 
by carrying out time integration in the half cycle on the crack 
opening side . 

Fig. 9 shows the relationship between AJc and evolution 
rate da/dN. Cracks are distributed along approximately the 
same linear line. Further, a multiplier term is near to 2 when 
it is arranged by the Paris law. It is contemplated that the 
initial crack behavior in the joint mode can be also evaluated 
using AJc. Fig. 9 also shows evolution characteristics (13) at 
room temperature that were arranged by the author using CCT 
specimen. It is supposed that a problem resides in the 
consistency of the steady evolution characteristics of the CCT 
specimen with that of bulk member, they have approximately the 
same tendency . 

da/dN = 32.1 xAJc x 807 ( 2 ) 
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Fig. 7 Analysis model 



Fig. 8 Creep J integrated value 

Fig. 9 Crack evolution behavior in solder bump joint 

4 . CONCLUSIONS 

The result of the mechanical fatigue test of solder bump 
joint, the prediction of crack initiation lifetime, and the 
simple prediction method of breakdown lifetime have been 
introduced. In addition to the above, the creep J integration 
range was determined by the FEM analysis, and the examined 
result of the crack evolution behavior in the actual joint mode 
has been reported, from which the following conclusions have 
been obtained. 

. In the lifetime prediction in the singularity region of 
non-elastic strain at bump angular point, a method of using the 
equivalent creep strain range, which is generated at a section 
apart from the singularity point by several grain sizes (50 pm) , 
as a parameter is applied. The result of lifetime prediction 
is in good correspondence to the result of test, and the 
reasonability of the method has been confirmed; 
. In the bump shape subjected to the test, the breakdown lifetime 
is about four to six times longer than the crack initiation 
lifetime. As the simple prediction method of bump breakdown 
lifetime in actual design, the authors have proposed the method 
that after the fatigue crack initiation lifetime is predicted 
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based on stress analysis, it is multiplied by the above lifetime 
ratio obtained by the test. 

. As the result that the crack evolution rate is evaluated by 
AJc calculated by the FEM analysis, there can be obtained a 
prospect that AJc can be arranged by the Paris law. There is 
also shown a possibility that the initial crack evolution 
behavior in actual joint mode is evaluated using AJc. 

. The above result was compared with the steady crack evolution 
behavior of bulk member obtained by CCT specimen. It is 
contemplated that there is a problem in the consistency between 
the steady evolution characteristics of bump joint and those 
of bulk member because the crack evolution behavior of bump 
joint is initial evolution behavior immediately after crack 
initiation at microscopic joint. However, in the scope of this 
study, both the bump joint and the bulk member have 
approximately the same tendency. 
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Appendix Calculation of creep J integrated value by FEM 
analysis - Basic examination using CCT specimen model - 

The creep J integrated value J* obtained by the FEM 
analysis was verified by a basic analysis model. The 
general-purpose structural analysis code ABAQUS was used. Fig. 
1 in Appendix shows the analysis model . A two-dimensional model 
in consideration of the symmetrical property of CCT specimen 
was formed by a plane stress element. The elastic-creep model 
was used as the solder material characteristics. Table 1 in 
Appendix shows creep characteristics expressed by a Norton law. 
Temperature conditions were set to constant temperatures of 25°C 
and -65°C. Forcible displacement 5 was applied to a load point 
in view of that many of solder joints were under displacement 
control type load conditions. Fig. 2 in Appendix shows J*s 
obtained in ten types of paths along the outer peripheries from 
the extreme end of a crack. Since J*s are initially in a 
small-scale creep state, J*s in the paths nearer to the extreme 
end of the crack are large. Thereafter, the specimen reaches 
a creep deformation managing type steady creep state in its 
entirety. J*s are converged to the same value although the 
accuracy of the J* in the first path is somewhat low, thereby 
independency of paths is exhibited. Fig. 2 in Appendix also 
shows J*s obtained by the following simple equation. 

n + 1 net [J) 
where a net shows net cross section nominal stress, and Vc shows 
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creep crack central opening displacement rate. In a small-scale 
creep state, the J*s obtained by the simple equation is 
particularly smaller than the value of the J* in the vicinity 
of the extreme end of the crack obtained by the FEM analysis. 
Further, under the condition of 25°C, the J*s obtained by the 
simple equation have a value about 16% smaller than that 
obtained by the FEM analysis. This is because a creep constant 
has a small value of n = 4.7, and thus the simple equation is 
insufficiently applied. Since n = 9.2 under the condition of 
-65°C, the difference between the J*s in the steady state is 
5% or less and they are approximately matched. 

Appended table Creep characteristics 

Appended Fig. 1 CCT specimen model 

Appended Fig. 2 Change of creep J integrated value 
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